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(54) Title: MULTI-DIMENSIONAL IMAGE RECONSTRUCTION 



(57) Abstract: Apparatus for radiation based imaging of a non-homogenous target area having distinguishable regions therein, com- 
prises: an imaging unit configured to obtain radiation intensity data from a target region in the spatial dimensions and at least one 
other dimension, and an image four-dimension analysis unit analyzes the intensity data in the spatial dimension and said at least one 
other dimension in order to map the distinguishable regions. The system typically detects rates of change over time in signals from 
radiopharmaceuticals and uses the rates of change to identify the tissues. In a preferred embodiment, two or more radiopharmaceu- 
ticals are used, the results of one being used as a constraint on the other. 
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Multi-Dimensional Image Reconstruction 



FIELD AND BACKGROUND OF THE INVENTION 
5 The present invention relates to multi-dimensional image reconstruction and, 

more particularly, but not exclusively to such image reconstruction based on a diffuse 
radioactive source or sources. 

Radiological imaging is generally carried out on a living target, which of 
course means a mix of tissues in close proximity, if not actually overlapping. The 

10 general procedure is to feed the patient with one or more radioactive markers prior to 
the imaging process. The radioactive markers are taken up by the digestive system 
and pass into the bloodstream. From the bloodstream the marker passes into the 
different tissues at varying rates depending on the tissue type. Some tissues absorb 
markers faster than others and some tissues absorb certain markers faster than others. 

15 Furthermore certain tissues flush out the markers faster than others, and again the rate 
of flushing out may also depend on the kind of marker being used. 

As a result, radioactive marking in fact creates a dynamic system in the body 
in which the relative darkness of a given tissue is related to a time factor. The 
radiologist knows that if he wants a good image of say the liver following application 

20 of a given marker then he should wait a certain number of hours from application of 
Ihe marker before taking the image. Even so, the liver is not dijfferentiated clearly 
from the other tissues. 

Examples of radiopharmaceuticals include monoclonal antibodies or other 
agents, e.g., fibrinogen or fluorodeoxyglucose, tagged with a radioactive isotope, e.g., 

25 99Mtechnetium, 67gallium, 201thallium, lllindiirai, 123iodine, 125iodine and 
ISfluorine, which may be administered orally or intravenously. The 
radiopharmaceuticals are designed to concentrate in the area of a tumor, and the 
uptake of such radiopharmaceuticals in the active part of a tumor, or other pathologies 
such as an inflammation, is higher and more rapid than in the tissue that neighbors the 

30 tumor. Thereafter, a radiation-emission-measuring-probe, which may be configured 
for extracorporeal or intracorporeal use, is employed for locating the position of the 
active area. Another application is the detection of blood clots with 
radiopharmaceuticals such as ACUTECT from Nycomed Amersham for the detection 
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of newly formed thrombosis in veins, or clots in arteries of the heart or brain, in an 
emergency or operating room. Yet other applications include radioimaging of 
myocardial infarct using agents such as radioactive anti-myosin antibodies, 
radioimaging specific cell types using radioactively tagged molecules (also known as 

5 molecular imaging), etc. 

The usual preferred emission for such applications is that of gamma rays, 
which emission is in the energy range of approximately 11-511 KeV, Beta radiation 
and positrons may also be detected. 

Radioactive-emission imaging is performed with a radioactive-emission- 

10 measuring detector, such as a room temperature, solid-state CdZnTe (CZT) detector, 
which is among the more promising that is currently available. It may be configured 
as a single-pixel or a multi-pixel detector, and may be obtained, for example, firom eV 
Products, a division of II- VI Corporation, Saxonburg Pa., 16056, or from IMARAD 
IMAGING SYSTEMS LTD., of Rehovot, ISRAEL, 76124, www.imarad.com. or 

15 from another source. Altematively, another solid-state detector such as CdTe, Hgl, 
Si, Ge, or the like, or a combination of a scintillation detector (such as Nal(Tl), LSO, 
GSO, Csl, CaF, or the like) and a photomultiplier, or another detector as known, may 
be used. 

Considering the issue in greater detail, certain biological or chemical 
20 substances such as targeted peptides, monoclonal antibodies and others, are used for 
tagging specific living molecules for diagnostic purposes. Ideally, these antibodies are 
specific to the desired type of cells, based on adhering only to specific molecular 
structures in which the antigene matching the antibody is highly expressed. The use of 
imaging devices such as a nuclear ganwna probe or a visual video probe can detect 
25 radiation emanating from taggants such as radionuclei or fluorescent dies that have 
been appended to the antibody before being delivered to the living body. An example 
is a cancerous cell of a prostate tumor on whose membrane there is an over expression 
of the Prostate Specific Membrane Antigen (PSMA). When a monoclonal antibody 
(Mab) such as Capromab Pendetide (commercially available as ProstaScint 
30 manufactured by Cytogen Corp.) is labeled with radioactive Indium (In 1 1 1) and is 
systemically delivered to the body, the Mab is carried by the blood stream and upon 
reaching the prostate tissue, adheres to the PSMA. The high energy radiation photons 
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emitted by the radioactive Indium can be detected using a nuclear camera, indicating 
the presence and the specific location of the tumor. 

Unfortunately, given the complexity of living organisms, in many instances 
the same antigen is also expressed in more than just the tissue under investigation. 
5 The antibody will thus also "paint" additional tissues such as infection areas, in 
addition to the tissue of interest. The radioactive readings taken from this additional 
tissue will be falsely interpreted as tumor areas, reducing the specificity of the test 
being performed. 

The 'Target to Backgroimd' ratio that characterizes every such antibody for a 
10 given target cell type is one of the major issues that determine the ability to perform 
proper diagnosis, and guided procedures. 

Since the uptake clearance of such a marker by the various tissues (target and 
background) varies over time, standard diagnosis protocols usually recommend taking 
an image at the time at which the ratio of Target emission vs. Background emission is 
15 the highest. 

In an experimental system tried out by researchers, two markers were supplied 
to various patients and then images were taken at successive intervals for each of the 
markers. Certain features in the target areas showed up clearly in all images, other 
features were clear for all images of one marker but faded in and faded out for the 
20 other marker, and yet other features faded in and out for both markers but at different 
times. The researchers were able to use their knowledge of the behaviors of the two 
markers with different tissues in order to identify the features in the images. 

The above system therefore relies on the knowledge of the researchers to put 
together information received from multiple images into an understanding of what the 
25 radio-imaging shows. In the general hospital environment it is not possible to 
guarantee that the necessary expertise is available, at least not for the amount of time 
that such a system would require. 

There is thus a widely recognized need for, and it would be highly 
advantageous to have, a radiological imaging system devoid of the above limitations. 

30 

SUMMARY OF THE INVENTION 
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According to one aspect of the present invention there is provided apparatus 
for radiation based imaging of a non-homogenous target area having distinguishable 
regions therein, the apparatus comprising: 

an imaging imit configured to obtain radiation intensity data from said target 
region in the spatial dimensions and at least one other dimension, and 

an image four-dimension analysis unit associated with said imaging unit for 
analyzing said obtained intensity data in said spatial dimension and said at least one 
other dimension in order to map said distinguishable regions. By image four- 
dimension analysis xmit is meant a unit that is able to analyze image data based on the 
spatial dimensions plus one other dimension, such as time. 

Preferably, said image four dimension analysis xmit is configured to constrain 
image output to a subset of said mapped distinguishable regions. 

Preferably, said image four-dimension analysis imit is configured to use said 
constraining to increase a resolution of said image output. 

Preferably, the distinguishable regions have differential takeup characteristics 
over time for radioactive markers and said multi-dimensional data includes a time 
component, said image analysis unit being configured to compare changes in detected 
intensity over time with said takeup characteristics in order to carry out said mapping. 

Preferably, at least two radioactive markers are applied to said target region, 
each of said markers having different takeup characteristics for respective ones of said 
regions, and each of said markers producing respectively distinguishable radiation, 
said image analysis unit being configured to use distinguishability of the radiation 
between said markers as an additional dunension in order to carry out said mapping. 

Preferably, at least two radioactive markers are applied to said target region, 
each of said markers having different takeup characteristics for respective ones of said 
regions, and each of said markers producing respectively distinguishable radiation, 
said image analysis unit being configured to use distinguishability of the radiation 
between said markers as said at least one other dimension in order to carry out said 
mapping. 

Preferably, said image analysis unit is configured to use said mapping to 
generate an image comprising said regions as distinct entities. 
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Preferably, said image analysis unit is configured to use said mapping to 
generate an image showing only a subset of said regions and to exclude at least one 
other of said regions. 

Preferably, said regions at least partially overlap, said image analysis unit 
being configured to show radiation due to a subset of said regions and to exclude 
radiation firom at least one other of said regions as noise. 

Preferably, said image analysis unit is configured to use said mapping to 
generate an image comprising said regions as distinct entities. 

Preferably, said image analysis unit is configured to use said mapping to 
generate an image showing only a subset of said regions and to exclude at least one 
other of said regions. 

Preferably, said regions at least partially overlap, said image analysis unit 
being configured to show radiation due to a subset of said regions and to exclude 
radiation fi-om at least one other of said regions as noise. 

Preferably, said imaging unit comprises at least one directional Geiger 
counter. 

Preferably, said imaging unit comprises a plurality of directional Geiger 
counters. 

Preferably, said imaging unit comprises a controller for directing said Geiger 
coimter to take images fi-om a set of locations optimized to obtain three-dimensional 
spatial data for a given target. 

Preferably, said non-homogenous target area is a region of living tissue, and 
said distinguishable regions are at least one of a group comprising: different tissues, 
different organs, blood and organ tissue, and tissue regions of differential pathologies. 

Preferably, one of said radioactive markers is thalium 201 and another of said 
radioactive markers is technetixmi 99. 

Preferably, said image analysis unit is configured to ignore image data as 
being firom outside said target area if said image data does not conform to at least one 
of said takeup characteristics. 

The apparatus may be configured to use said mapping to identify at least one 
region of low emissivity, thereby to concentrate imaging resources on said identified 
region. 
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Preferably, said imaging is via voxels of s^d target area and wherein said 
concentrating imaging resources comprises merging voxels of said identified region. 

Preferably, said concentrating resources comprises concentrating resources of 
said imaging unit on said identified region. 

Preferably, said mapping comprises a first mapping to identify an organ and a 
second mapping constrained within said organ. 

According to a second aspect of the present invention, there is provided a 
method for constraining image data obtained fi-om a non-homogenous target area 
having a plurality of distinct regions using at least one radioactive marker, the method 
comprising: 

obtaining radiation intensity data fi:om said target area m the spatial 
dimensions and at least one other dimension, and 

analyzing said radiation intensity data using said dimensions in order to 
classify data pertaining to respective distinct regions and to produce an output 
constrained to a subset of said distinct regions. 

Preferably, said at least one radioactive marker has a time absorption 
characteristic which is different for ones of said distinct regions and wherein said at 
least one other dimension is a time dimension. 

Preferably, there are provided at least two radioactive markers, each having 
respectively different time absorption characteristics for each of said distinct regions, 
and wherein said at least one other dimension is a time dimension. 

Preferably, there are provided at least two radioactive markers, each producing 
a distinct signal and each respectively having different time absorption characteristics 
for each of said distinct regions, and wherein said at least one other dimension is a 
dimension indicating individual markers. 

Preferably, there are provided at least two radioactive markers, each producing 
a distinct signal and each respectively having different time absorption characteristics 
for each of said distinct regions, and wherein said at least one other dimension 
comprises a time dimension and a dimension indicating individual markers. 

Preferably, said classifying identifies at least one region of low emissivity, the 
method further comprising concentrating imaging resources on said identified region. 

Preferably, said imaging is via voxels of said target area and said 
concentrating imaging resources comprises merging voxels of said identified region. 
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Preferably, said concentrating resources comprises concentrating resoxirces of 
said imaging unit on said identified region. 

Preferably, said classifying comprises a first classification to identify an organ 
and a second classification constrained within said organ. 

According to a third aspect of the present invention there is provided apparatus 
for radiation based imaging of a non-homogenous target area having distinguishable 
regions therein, the apparatus comprising: 

an imaging unit configured to obtain radiation intensity data from said target 
region in the spatial dimensions and at least one other dimension, 

an image analysis unit associated with said imaging unit for analyzing said 
obtained intensity data in said spatial dimension and said at least one other dimension 
in order to map said distinguishable regions, and using said mapping to control use of 
imaging resources. 

Unless otherwise defined, all technical and scientific terms used herein have 
the same meaning as commonly understood by one of ordinary skill in the art to 
which this invention belongs. The materials, methods, and examples provided herein 
are illustrative only and not intended to be limiting. 

Implementation of the method and system of the present invention involves 
performing or completing certain selected tasks or steps manually, automatically, or a 
combination thereof. Moreover, according to actual instrumentation and equipment 
of preferred embodiments of the method and system of the present invention, several 
5 selected steps could be implemented by hardware or by software on any operating 
system of any firmware or a combination thereof- For example, as hardware, selected 
steps of the invention could be implemented as a chip or a circuit. As software, 
selected steps of the invention could be implemented as a plurality of software 
instructions being executed by a computer using any suitable operating system. In 
10 any case, selected steps of the method and system of the invention could be described 
as being performed by a data processor, such as a computing platform for executing a 
plurality of instructions. 

BRIEF DESCRIPTION OF THE DRAWINGS 
15 The invention is herein described, by way of example only, with reference to 

the accompanying drawings. With specific reference now to the drawings in detail, it 
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is stressed that the particulars shown are by way of example and for purposes of 
illustrative discussion of the preferred embodiments of the present invention only, and 
are presented in order to provide what is believed to be the most useful and readily 
xmderstood description of the principles and conceptual aspects of the invention. In 
5 this regard, no attempt is made to show structural details of the invention in more 
detail than is necessary for a fundamental understanding of the invention, the 
description taken with the drawings making apparent to those skilled in the art how 
the several forms of the invention may be embodied in practice. 
In the drawings: 

10 FIG, 1 is a simplified diagram showing a single detector detecting over a 

target region; 

FIG. 2 is a simplified diagram showing two detector positions (not necessarily 
simultaneously) allowing three-dimensional information to be obtained from a target 
region; 

15 FIGS. 3A-3D show a series of four time absorption characteristics for different 

radiopharmaceuticals within different tissues; 

FIG. 4 is a simplified schematic diagram showing a device for driving an 
imaging head and allowing control of the imaging head by the image analyzer device; 
FIG. 5 is a simplified flow chart illustrating the image analysis process carried 
20 out by the analyzer in Fig. 4 in the case of a single marker; 

FIGS. 6A-6D illustrate two sets of successive images of the same target area 
taken using two different markers respectively, according to a preferred embodiment 
of the present invention; 

FIG. 7A is a simplified flow chart illustrating a procedure according to a 
25 preferred embodiment of the present invention using two or more markers for first of 
all identifying an organ and then secondly determining the presence or otherwise of a 
pathology within that organ; 

FIG. 7B is a simplified flow chart showing a generalization of Fig 7 A for the 
general case of two specific pattems; 
30 FIG. 8 is a simplified flow chart illustrating a procedure according to a 

preferred embodiment of the present invention using two or more markers for 
identifying a region of low emissivity within a target area and using that identification 
to control imaging resources to better image the identified region; and 
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FIGS. 9 A - 9D illustrate two sets of successive images of the same target area 
taken using two different markers, in a similar way to that shown in FIG. 6, except 
that this time the regions of interest are one inside the other. 

5 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present embodiments comprise an apparatus and a method for radiation 
based imaging of a non-homogenous target area having regions of different material 
or tissue type or pathology. The imaging uses multi-dimensional data of the target 
area in order to distinguish the different regions. Typically the multi-dimensional 

10 data involves time as one of the dimensions. A radioactive marker has particular 
time-absorption characteristics which are specific for the different tissues, and the 
imaging device is programmed to constrain its imaging to a particular characteristic. 

The result is not merely an image which concentrates on the tissue of interest 
but also, because it is constrained to the tissue of interest, is able to concentrate 

15 imaging resources on that tissue and thus produce a higher resolution image than the 
prior art systems which are completely tissue blind. 

The principles and operation of a radiological imaging system according to the 
present invention may be better imderstood with reference to the drawings and 
accompanying description. 

20 Beifore explaining at least one embodiment of the invention in detail, it is to be 

understood that the invention is not limited in its application to the details of 
construction and the arrangement of the components set forth in the following 
description or illustrated in the drawings. The invention is capable of other 
embodiments or of being practiced or carried out in various ways. Also, it is to be 

25 vmderstood that the phraseology and temiinology employed herein is for the pvirpose 
of description and should not be regarded as limiting. 

Reference is now made to Fig. 1, which illustrates a simple Geiger coimter 
taking an image of a target according to the prior art. Geiger counter 10 is placed in 
association with target 12 and absorbs any radioactive particles that come its way. In 
general the radioactive particles arriving at the Geiger counter arrive from somewhere 
within cone 14. The Geiger coxmter has no information as to the depth from which 
the particle comes and cannot even distinguish between particles arriving from 
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different directions within the cone. Thus in principle the prior art Geiger counter 
gives low resolution one dimensional information. 

If the counter is now moved to different positions over the surface of the target 
then the data from the different positions can be built up into a low resolution two- 
dimensional image. 

One way of increasing the resolution of the Geiger counter is to make it 
smaller. Then the cone, whilst retaining the same geometry, gives higher resolution 
data. 

The detector takes (yt)^^i samples to form a data set, which would typically be 

a two-dimensional image of the target from a given direction. 

Reference is now made to Fig. 2, which is a simplified diagram showing how 
three-dimensional information can be obtained from the target. Parts that are the same 
as in previous figures are given the same reference numerals and are not referred to 
again except as necessary for understanding the present embodiment. A second 
Geiger counter 16 is placed essentially at right angles to the first Geiger coxmter and 
obtains a similar kind of image to the first Geiger coxmter. However, since the two 
cones overlap, the images produced can be cross-correlated to infer the presence of 
hot or cold radiation sources in three dimensions. 

Reference is now made to Fig. 3, which is a sequence of graphs illustrating the 
different absorption characteristics for different tissues of a given radioactive marker. 
Typical markers that may be considered are Thalium 201 and Technitium 99. Fig. 3a 
indicates a typical absorption characteristic of thalium 201 for blood, thalium 201 
being a particularly good marker for blood. The marker is generally absorbed by the 
blood fairly rapidly following digestion and then gradually disappears as it is taken up 
by the various tissues and organs including the kidneys. Marker material from the 
tissues eventually finds its way back into the blood for excretion. That which is 
absorbed by the kidneys is excreted directly and not seen again. 

Figures 3B, 3C and 3D show time absorption characteristics for technitium 99 
for different tissues, and it will be seen that the characteristic is generally curved but 
peaks at different times for the different tissues. 

The principle on which the present embodiments are based is as follows: 
Considering the graphs in Fig. 3, it will be apparent that a region belonging to a single 
tissue will behave in a uniform manner as regards signal intensity. That is to say, a 
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given marker will be taken up and then expelled at the same rate over a given tissue, 
whereas this rate will be different for other tissues. If therefore a series of successive 
images are taken of the target and the images are analyzed region by region for rates 
of change of intensity, a particular desired region can be identified by virtue of having 
rates of change in intensity that fit with a given characteristic. The regions are 
distinguishable in this way even if the region of interest is heavily overlapped with 
other regions. 

Reference is now made to Fig. 4, which shows apparatus for radiation-based 
imaging of a non-homogenous target area. Apparatus 20 comprises an imaging unit 
22 which itself consists of a series of small Geiger counters 24.1...24.n arranged on 
an imaging head. The imaging unit is controlled by motion controller 26 to take 
readings from different locations around the target area. Preferably, the motion of the 
imaging head is controlled by software via servo-motors. In addition the motions, 
either of the individual Geiger counters or of groupings of the Geiger coxmters, is also 
controlled by software via servo-motors. 

In a preferred embodiment, the signals received from the individual Geiger 
counters are summed to fomi a three-dimensional image of the target area. The 
skilled person will appreciate that the system could also be based on a two- 
dimensional image. In either case, the signals are fed to an image analyzer 28, where 
the signals are analyzed to form images. 

In the preferred embodiments, the image analyzer is able to use the marker 
take up characteristics to compare successive images and identify regions of particular 
interest, and then to concentrate imaging resources on tiiose regions. That is to say 
the image analyzer is in fact able to control further operation of the imager. 

Reference is now made to Fig. 5, which is a simplified flow chart illustrating 
the image analysis process carried out by analyzer 28 in the case of a single marker. 
Preferably a series of images of the same views are taken at different times, stage 30, 
and a three-dimensional overall image of the target is formed for each time. The 
analyzer then analyzes each of the three-dimensional overall images for local 
intensities at different locations aroxmd the target, stage 32. The local intensities are 
noted and the same locations on the different images are superimposed in stage 34. 
From the superpositioning, local rates of change of intensity between the images may 
be obtained in stage 36. The rates of change are compared with the pre-obtained 
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characteristics for the marker with the different tissues in stage 38, and the data are 
then constrained to those localities which conform to the desired predetermined 
characteristics in stage 40. As a result the imaging process can be used to identify and 
concentrate on localities of interest and data from other localities can be jettisoned. 
Consequently, the image analysis is able to concentrate its resources on the tissues of 
interest and a higher resolution final image can be produced. 

It will be appreciated that in many cases two types of tissue may be 
superimposed, of which only one of the tissues is of interest In this case it is of equal 
importance both to exclude the one tissue that is not of interest and to include the 
tissue that is of interest It may be that the best marker for one tissue may not be the 
best marker for the other tissue. The system as described with respect to Figs 4 and 5 
may be adapted to use with two or more markers, as exemplified in Fig. 6. Each 
marker produces a radioactive particle of different energy level, and therefore the data 
from the different markers can be collected and summed separately to form different 
images. Mathematically the different data sets obtained from the different energy 
level signals may be treated as different dimensions of a multi-dimensional vector. 
For each of the marker-images the appropriate characteristics are used to identify the 
tissues of interest, and the results can be cross-checked between the different markers. 
The different tissues can be mapped and the image analysis can concentrate on the 
area of interest As a result the system uses both time and particle energy as separate 
dimensions in addition to the spatial dimensions in order to characterize or map the 
tissues. 

As a result the image analysis unit is able to produce a final result treating the 
various tissue regions as separate entities. Furthermore, as the system is aware of the 
regions as entities it is able to further direct the imaging process to concentrate on the 
regions of interest. 

An example in which regions at least partially overlap is the heart. Generally, 
scans of the heart are interested in the muscular walls of the heart. Although the 
chambers of the heart are filled with blood, any signal coming from the blood is in 
fact noise to this kind of scan. It is therefore advantageous to carry out an imaging 
process which is able to positively identify signals firom the muscular heart walls and 
at the same time exclude the blood. 



wo 2005/067383 



PCT/IL2005/000048 



13 

Referring now to Fig. 6, and in a preferred embodiment, the patient ingests 
two markers, thalium 201 and technetium 99. The first of these is an effective blood 
marker and two successive thalium images are shown in Figs 6a and 6b, and the 
second is more effective at marking muscle tissue and two successive images thereof 
are shown in Figs 6c and 6d. The heart is imaged at intervals chosen both for the 
characteristic for thalium 201 in blood and for the characteristic of technetium 99 in 
muscle. The result is a series of images for each of the markers. The series for 
thalium 201 may be constrained to show the regions of blood quite clearly, and to 
filter out other regions. In here a blood vessel is shown clearly in 6a and more faintly 
in 6b where the thalium has mostly been flushed out. The series for technetium 99, 
figs 6c and 6d show muscle wall structures. The first of the two images apparently 
shows larger structures but in fact all that it is showing is that much technetium has 
not yet been absorbed in the muscle. The second image 6d may therefore be used to 
constrain the first image 6c to show only the muscle walls regions. The two series of 
images may then be superimposed to filter out from the technetium 99 images 6c and 
6d anj^hing that appears strongly in the thaUimi images 6a and 6b, The filtering may 
additionally remove anything that appears strongly in both images as coming fiom 
outside the region. 

In the above example, two regions were of respectively positive and negative 
interest, meaning one for concentrating on and the other for filtering out. It will be 
appreciated that several regions or several tissue types may be of positive interest or 
there may be any combination of regions with just one being of positive interest. 
Alternatively all regions may be of positive interest but importance may be attached 
to discriminating between the different signals from the different regions. 

The system is able to use the mapping to generate an image comprising the 
different tissue regions as distinct entities. As a consequence of the mapping process, 
the system is able to be aware electronically of the different regions and thus control 
both the imaging head and the analysis unit to concentrate their resources on specific 
regions. The result is greater resolution for the regions of interest. 

The preferred embodiments may be used to expand the information obtained 
from the markers, using either or both of examining the kinetics of the markers over 
time and using several markers concurrently. 
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In order to increase the specificity of the test, additional second substances 
("secondary substances"), with reactivity and pharmaco-kinetics differing from those 
of the first substance can be used in order to enhance the differentiation between the 
different pathologies, as explained above with respect to Fig. 6, The secondary 
substance, in this case thalium, ideally marks only a subset of the population marked 
by the primary substance and does so at different rates. Such a difference exists 
because of different affinity to various cell types and different participation in 
metabolic reactions of different tissues. The difference is associated with the rate of 
marking and/or with the location of the marking. 

Upon reading the radioactive signals emanating from the voxels stemming 
from different substances at different time instances, it is possible to build for every 
voxel a multi dimensional data matrix Sjk whose elements are intensity readings taken 
at instances K resulting from the interaction of Substance J. Examination of every 
voxel of tissue in this multidimensional space quantifies the temporal and specific 
reaction of the tissue to different substances and thus increases the probability of 
specific detection of different pathologies. Furthermore, standard image processing 
techniques can be used in order to more acciirately define the spatial location of 
different pathologies. 

In addition to the method above, spatial properties that reflect typical 
relationships between neighboring voxels may also be a criteria and represented as 
part of the pattern of the tissue type. 

Reference is now made to Fig. 7, which illustrates an additional statistical 
approach. In Fig. 7, an automatic algorithm based on expected intensities may be 
used to determine if the entire organ or region is diseased or non-diseased. Once it is 
possible to become tissue-aware, as explained above, then it is no longer necessary to 
carry out such analysis on a voxel-by- voxel basis. Rather the system is able to 
determine where the organ lies say using a first marker and then a second marker may 
be imaged using the constraint of the organ location, the second marker being able to 
locate the presence of the pathology. 

Reference is now made to Fig. 8 which illiistrates a method for using the tissue 
aware properties of the present embodiments in order to tune detection to match tissue 
or organ emissivities. Generally, any region, no matter how much radiation it 
produces, can always be imaged sufficiently simply by leaving the measuring device 
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in position for long enough. However, in many cases there may be limited time 
available. For such cases in which there is limited time for data acquisition, the 
present embodiments can be used to identify regions that may be expected to produce 
less emission. The system may then tune imaging resources or resolution onto those 
tissues according to the nxraiber of photons available. Clearly the more photons 
obtained the more reliable is the data, and therefore a tissue aware system is able to 
concentrate more detectors on the weaker signaling tissues. 

If there are still not enough photons, or there are not enough detectors, then 
another way of pooling resources is to merge neighboring voxels (or regions). Such a 
procedure may reduce resolution, but will increase the overall number of photons for 
that merged region, and thus enable better classification of that region based on a 
more reliable photon count. Such a compromise enables analysis of the same 
collected data by ways that would allow high resolution where there are enough 
photons and lower resolutions where there are less while maintaining reliability of the 
analysis. 

Again the tissue regions may be identified using multiple markers. 

The above-described embodiment may lead to controlled sensitivity levels, 
currently not available with radioimaging. 

The concept of using multiple antibodies can be used for therapy purposes, as 
in the following: 

The specificity of a single antibody carrying a drug (or radioactive 
therapy) determines the chance for non-target tissue to receive the drug, and thus be 
subject to any toxicity of the drug. In cases where there are several antibodies, each 
with limited specificity, but with affinity to different 'background' tissue, a 
combination of antibodies may be used to improve the overall specificity, and thus to 
reduce overall toxicity and enable higher efficacy of treatment. 

For example, if a first antibody (Al) based dmg binds to the target Nl folds its 
affinity to the closest non-target tissue (Bl), and a second antibody (A2) with similar 
dmg has target affinity which is N2 folds higher than its closest non-target tissue 
(B2), then using a merged therapy will enable better target vs. non-target specificity, 
which is better than Nl and N2 (assuming Bl and B2 are different). 

In a more generalized embodiment, the system may include a signal analysis 
module, including a library of patterns that are typical for every cell type. Each type 
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of cells has one or more patterns associated with it, and the pattern determines how a 
set of markers injected according to a specific protocol (dosage, time, etc) may be 
expressed in that cell type. The analysis includes classification of the readings from 
each voxel based on correlation, or other statistical tools for assessing the most 
probable tissue classification for each voxel. 

Since there may be several cell types for a given disease (e.g. cancer may 
show in several forms), the algorithm may be optimized to determine the exact tissue 
type per voxel or region. Altematively, the algorithm may be optimized to determine 
the general property of diseased/non-diseased, while taking the specific classification 
only as a factor in the statistical analysis. 

It should be noted that the system may allow for various protocols for 
administering the markers, where injection of the various markers may be 
simultaneous, or multiple injections at various times, as various markers have 
different lifetime in the circulation. 

The issue of generating imaging using two or more markers is now considered 
mathematically. 

An intensity distribution /, conventionally defined in terms of radioactive 
emissions per seconds, is now redefined as a vector of distributions over the volume 
U, forming our input space. Each dimension of the vector is a different one of the 
radiopharmaceuticals. The universal set U comprises a set of basic elements u {e.g., 
pixels in two dimensional spaces, voxels in three dimensional spaces), and I{u) is the 
intensity in a given basic element u g U. For j radiopharmaceuticals this becomes 
I(uf^^^ An inverse (or reconstruction) problem arises when one caimot sample directly 
from /, but can sample from a given set of views <S>. A projection (pe€> is defined by 
the set of probabilities {(p (u):ugU}, where cp (u) is the probability of detecting a 
radioactive emission from a voxel u, as defined by viewing parameters, such as the 
physical and geometrical properties of the detecting unit, as well as the attenuation 
parameters of the viewed volume U, and the time parameters. A measurement is 
obtained by choosing a view <peO, and then sampling according to the viewing 
parameters. 

For j radiopharmaceuticals or markers and k detectors, the probability of 
seeing a particle becomes (pj^ (w) 
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In the following analysis, / is the intensity of a radioactive substance, and the 
viewing parameters include the geometrical properties of a coUimated detecting unit 
and the detecting unit's position and orientation with respect to volume U. The 
number of radioactive emissions counted by the detecting unit within a time interval 
is a Poisson distribution, where <p (u) is the detection probability of a photon emitted 
from voxel u sU and the mean of the distribution is the weighted sum S„^£/ cp (u)I(u), 

For the case of the kth detector a measurement Yk = ll>ueuXt(u)^ where X(U) is 
a Poisson distribution. 

Where Y^xt) = ^^OM (m). 
Hence Y(j,k,t)) = Poisson (Y(j,k,t)) 

The projection set is thus defined by a matrix whose rows are the 
projections of the chosen views. 7 is a vector of densities (specified per each element 
in U)^ and O/ is a vector of respective effective intensity levels for the views in the 
set. A vector of measurements y is obtained by a random sample from each view 
(according to the associated Poisson distribution). As discussed above, there are 
various known reconstruction methods that provide estimators for / given the 
projections O and the measurements y. 

Using the above mathematics the problem is solved (an image created) one of 
the vectors say once an hour. The rates of change are determined. Simultaneously 
the problem is solved for another of the vectors at similar time intervals and the rates 
of change are determined. Then a stage of cross-identijScation is carried out between 
the two images, so that wanted tissues as identified by each image minus imwanted 
tissues identified by each image are concentrated on to form a new image. Cross- 
identification may be an iterative process. 

In the example given above of the imaging of the heart using one blood 
marker and one muscular tissue marker, the areas identified by the blood marker are 
subtracted. The areas identified by the muscle marker are added, and those tissues not 
identified by either are likewise ignored as being signals from outside the target 
region. 

The non-homogenous target area is typically a region of living tissue, 
generally belonging to a patient. The distinguishable regions within can be different 
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tissues, different organs, a mixture of blood and organ tissue as witii the above 
example of the heart, or tissue regions exhibiting differential pathologies. 

It is expected that during the life of this patent many relevant markers, 
radiological imaging devices and two and three dimensional imaging systems will be 
developed and the scopes of the corresponding terms herein, are intended to include 
all such new technologies a priori. 

5 It is appreciated that certain features of the invention, which are, for clarity, 

described in the context of separate embodiments, may also be provided in 
combination in a single embodiment. Conversely, various features of the invention, 
which are, for brevity, described in the context of a single embodiment, may also be 
provided separately or in any suitable subcombination. 

10 Although the invention has been described in conjunction with specific 

embodiments thereof, it is evident that many alternatives, modifications and variations 
will be apparent to those skilled in the art. Accordingly, it is intended to embrace all 
such altematives, modifications and variations that fall within the spirit and broad 
scope of the appended claims. All publications, patents and patent applications 

15 mentioned in this specification are herein incorporated in their entirety by reference 
into the specification, to the same extent as if each individual publication, patent or 
patent application was specifically and individually indicated to be incorporated 
herein by reference. In addition, citation or identification of any reference in this 
application shall not be construed as an admission that such reference is available as 

20 prior art to the present invention. 



wo 2005/067383 



PCT/IL2005/000048 



19 

WHAT IS CLAIMED IS: 

1. Apparatus for radiation based imaging of a non-homogenous target 
area having distinguishable regions therein, the apparatus comprising: 

an imaging unit configured to obtain radiation intensity data from said target 
region in the spatisd dimensions and at least one other dimension, and 

an image four-dimension analysis unit associated with said imaging unit for 
analj^ing said obtained intensity data in said spatial dimension and said at least one 
other dimension in order to map said distinguishable regions. 

2. Apparatus according to claim 1, wherein s^d image four dimension 
analysis unit is configured to constrain image output to a subset of said mapped 
distinguishable regions. 

3. Apparatus according to claim 2, wherein said image four-dimension 
analysis unit is configured to use said constraining to increase a resolution of said 
image output. 

4. Apparatus according to claim 1, wherein said distinguishable regions 
have differential takeup characteristics over time for radioactive markers and said 
multi-dimensional data includes a time component, said image analysis unit being 
configured to compare changes in detected intensity over time with said takeup 
characteristics in order to carry out said mapping. 

5. Apparatus according to claim 4, wherein at least two radioactive 
markers are applied to said target region, each of said markers having different takeup 
characteristics for respective ones of said regions, and each of said markers producing 
respectively distinguishable radiation, said image analysis unit being configured to 
use distinguishability of the radiation between said markers as an additional 
dimension in order to carry out said mapping. 

6. Apparatus according to claim 1, wherein at least two radioactive 
markers are applied to said target region, each of said markers having different takeup 
characteristics for respective ones of said regions, and each of said markers producing 
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respectively distinguishable radiation, said image analysis unit being configured to 
use distinguishability of the radiation between said markers as said at least one other 
dimension in order to carry out said mapping. 

7. Apparatus according to claim 4, wherein said image analysis imit is 
configured to use said mapping to generate an image comprising said regions as 
distinct entities. 

8. Apparatus according to Claim 4, wherein said image analysis unit is 
configured to use said mapping to generate an image showing only a subset of said 
regions and to exclude at least one other of said regions. 

9. Apparatus according to claim 4, wherein said regions at least partially 
overlap, said image analysis unit being configured to show radiation due to a subset of 
said regions and to exclude radiation from at least one other of said regions as noise. 

10. Apparatus according to claim 5, wherein said image analysis unit is 
configured to use said mapping to generate an image comprising said regions as 
distinct entities. 

11. Apparatus according to claim 5, wherein said image analysis xmit is 
configured to use said mapping to generate an image showing only a subset of said 
regions and to exclude at least one other of said regions. 

12. Apparatus according to claim 5, wherein said regions at least partially 
overlap, said image analysis unit being configured to show radiation due to a subset of 
said regions and to exclude radiation from at least one other of said regions as noise, 

13. Apparatus according to claim 1, wherein said imaging xmit comprises 
at least one directional Geiger counter, 

14. Apparatus according to claim 13, wherein said imaging unit comprises 
a plurality of directional Geiger counters. 
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15. Apparatus according to claim 13, wherein said imaging unit comprises 
a controller for directing said Geiger counter to take images from a set of locations 
optimized to obtain three-dimensional spatial data for a given target, 

16. Apparatus according to claim 1, wherein said non-homogenous target 
area is a region of living tissue, and said distinguishable regions are at least one of a 
group comprising: different tissues, different organs, blood and organ tissue, and 
tissue regions of differential pathologies. 

17. Apparatus according to claim 5, wherein one of said radioactive 
markers is thalium 201 and another of said radioactive markers is technetium 99. 

18. Apparatus according to claim 4, wherein said image analysis unit is 
configured to ignore image data as being from outside said target area if said image 
data does not conform to at least one of said takeup characteristics. 

19. Apparatus according to claim 1, configured to use said mapping to 
identify at least one region of low emissivity, thereby to concentrate imaging 
resources on said identified region, 

20. Apparatus according to claim 19, wherein said imaging is via voxels of 
said target area and wherein said concentrating imaging resoiurces comprises merging 
voxels of said identified region. 

21. Apparatus according to claim 19, wherein said concentrating resources 
comprises concentrating resources of said imaging unit on said identified region. 

22. Apparatus according to claim 1, wherein said mapping comprises a 
first mapping to identify an organ and a second mapping constrained within said 
organ. 
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23. A method for constraining image data obtained from a non- 
homogenous target area having a piuraUty of distinct regions using at least one 
radioactive marker, the method comprising: 

obtaining radiation intensity data from said target area in the spatial 
dimensions and at least one other dimension, and 

analyzing said radiation intensity data using said dimensions in order to 
classify data pertaining to respective distinct regions and to produce an output 
constrained to a subset of said distinct regions. 

24. The method of claim 23, wherein said at least one radioactive marker 
has a time absorption characteristic which is different for ones of said distinct regions 
and wherein said at least one other dimension is a time dimension. 

25. The method of claim 24, wherein there are provided at least two 
radioactive markers, each having respectively different time absorption characteristics 
for each of said distinct regions, and wherein said at least one other dimension is a 
time dimension. 

26. The method of claim 23, wherein there are provided at least two 
radioactive markers, each producing a distinct signal and each respectively having 
different time absorption characteristics for each of said distinct regions, and wherein 
said at least one other dimension is a dimension indicating individual markers. 

27. The method of claim 23, wherein there are provided at least two 
radioactive markers, each producing a distinct signal and each respectively having 
different time absorption characteristics for each of said distinct regions, and wherein 
said at least one other dimension comprises a time dimension and a dimension 
indicating individual markers. 

28. The method of claim 23, wherein said classifying identifies at least one 
region of low emissivily, the method further comprising concentrating imaging 
resources on said identified region. 
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29. The method of claim 28, wherein said imaging is via voxels of said 
target area and wherein said concentrating imaging resources comprises merging 
voxels of said identified region. 

30. The method of claim 28, wherein said concentrating resources 
comprises concentrating resources of said imaging unit on said identified region. 

31. The method of claim 23, wherein said classifying comprises a first 
classification to identify an organ and a second classification constrained within said 
organ. 

32. Apparatus for radiation based imaging of a non-homogenous target 
area having distinguishable regions therein, the apparatus comprising: 

an imaging unit configured to obtain radiation intensity data from said target 
region in the spatial dimensions and at least one other dimension, 

an image analysis unit associated with said imaging unit for analyzing said 
obtained intensity data in said spatial dimension and said at least one other dimension 
in order to map said distinguishable regions, and using said mapping to control use of 
imaging resources. 
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